Decoupling of Spurious Deep Bound States 
with the Similarity Renormalization Group 



K.A. WendtQ R.J. FurnstahlQ and R.J. Perrj0 

Department of Physics, The Ohio State University, Columbus, OH 43210 
(Dated: January 17, 2011) 

The Similarity Renormalization Group (SRG) is a continuous series of unitary transformations 
that can be implemented as a flow equation. When the relative kinetic energy (Ti-ci) is used in the 
SRG generator, nuclear structure calculations have shown greatly improved convergence with basis 
size because of the decoupling of high-energy and low-energy physics. However this generator can 
sometimes be problematic. A test case is provided by a study of initial interactions from chiral 
effective field theories with large cutoffs, which can lead to spurious deep bound states. We would 
like the SRG to decouple these from the physical shallow bound states. However, with Troi the 
high- and low-energy bound states are not decoupled in the usual sense. Replacing Trci by the 
momentum-space diagonal of the Hamiltonian (Hd) in the SRG generator does produce decoupling, 
such that the shallow states are in the low-momentum region and the deep bound states are at 
higher momentum. The flow toward universal low-momentum interactions is also restored. 



I. INTRODUCTION 



where 



Similarity Renormalization Group (SRG) flow equa- 
tions with a generator based on the relative kinetic en- 
ergy have been shown to soften intcr-nucleon interactions 
by decoupling energy scales [U [5] , leading to significantly 
enhanced convergence in calculations of nuclear structure 
and reactions [2HS]- To date, these SRG applications 
have all started with phenomenological potentials such 
as Argonne vis [5] or the N^LO chiral potentials from 
Refs. [10] and [H]. When evolved to softened form, ma- 
trix elements of these interactions are found to collapse 
toward a common NN potential for momenta below the 
decoupling scale [2]. These desirable SRG features are 
not guaranteed, however, when applied to chiral effective 
field theory (EFT) interactions with large cutoffs. Such 
interactions have been explored at leading order [T2] to 
examine renormalization issues (see also Refs. [T3HI1]). 
The SRG approach as applied so far to nuclear interac- 
tions fails for these theories in channels where the tensor 
force from pion exchange introduces spurious deep bound 
states. The generic cause of this failure was identified 
by Glazek and Perry [20]. In this paper we document 
the specific problems with large-cutoff chiral EFT (xeft) 
and test whether they are fixed as described in Ref. [5D] . 

The SRG as used here is implemented as a series of 
infinitesimal unitary transforms of a Hamiltonian H, la- 
beled by a fiow parameter s [5TJ [52] , 



H{s) = U{s)H{s = 0)U'f{s) 



(1) 



Differentiating with respect to the flow parameter shows 
how any evolution is constrained by unitarity. 



ds 



H{s)^[fjis),H{s)] , 



(2) 
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ry(s) 



dUis] 



(3) 



The fiow is controlled through the choice of the anti- 
Hermitian generator ri(s). To generate a renormalization 
group evolution, we need the flow to decouple high and 
low energy degrees of freedom (i.e. to remove far off diag- 
onal matrix elements). For many applications, a useful 
form for the generator is 



[G{s),His)] 



(4) 



where G{s) is a Hermitian operator. The original choice 
for G{s) advocated by Wegner and collaborators |2T1 [22] 
and applied extensively in condensed matter is the diag- 
onal component of the interaction, G{s) — Hd{s), 



if « = i , 
otherwise. 



(5) 



If the basis in which is taken to be diagonal is a 
discretized, partial- wave momentum basis \i) = |fcjm), 
the Hamiltonian is driven toward diagonal form with s 
such that the parameter A = s^^/^ is a measure of the 
width of diagonal dH [52] . 

In most applications to nuclear systems, G{s) has been 
taken simply to be the relative kinetic energy T^oi, which 
is by construction independent of s. This choice also has 
the effect in practice of driving the momentum-space po- 
tential toward the diagonal, which directly softens the 
repulsive core of nucleon-nucleon interactions (as well 
as short-range tensor forces) by decoupling high- and 
low- momentum degrees of freedom [l]. No adverse re- 
sults have been observed for light nuclei; however, this 
SRG evolution has always been terminated such that 
A > 1 fm-\ 

In Ref. [2^ , it was observed that when evolving a sim- 
ple model Hamiltonian, the Wegner evolution (G = Hd) 
will decouple the bound state by leaving it as an isolated 



2 



eigenvalue on the diagonal of the Hamiltonian matrix. 
In contrast, they found a very different behavior with 
the evolution generated using T^^i- With Treb the bound 
state remains coupled to low momentum, and is pushed 
to the lowest momentum part of the matrix. This unde- 
sirable result has not been seen in the two-nucleon system 
for A > 1 fm~^ because the only physical bound state, 
the deuteron, is very shallow (e.g., the scattering length 
in this channel is unnaturally large) and only plays a role 
for A <C 1 fm~^. Indeed, the evolution for T^ci and Hd in 
the usual range of A is practically identical numerically 
(at least for the A = 2 system). 

In most cases, we do not want to decrease A to the 
point where bound-state eigenvalues appear on the diag- 
onal. The effects of discretized momenta are magnified 
when H is driven too close to diagonal form, so we typi- 
cally seek a window in which high-energy states decouple 
but the low-energy part of H is insensitive to discretiza- 
tion. However, if a deep spurious bound state exists be- 
cause our effective Hamiltonian is singular at short dis- 
tances, we want to make sure that any associated spu- 
rious physics decouples from the low-energy physics our 
effective theory is intended to capture. In this way the 
practical advantages of the SRG are preserved. This is 
accomplished if the spurious bound eigenstate is forced to 
the diagonal when A is still sufficiently above the scale of 
interest. According to Glazek and Perry, this happens for 
the iJjj-based generator but T^^i drives the bound-state 
eigenvalue to the low-energy corner of the Hamiltonian 
matrix |20) . As a result, decouples the spurious state 
from low-energy physics but Tj-ei causes the spurious state 
to corrupt low-energy physics. 

In this paper, we test whether the analysis of Glazek 
and Perry carries over to the nuclear case with spuri- 
ous bound states and further investigate the nature of 
decoupling and the flow toward universal interactions in 
the SRG. In Section [Hj we give background details on 
the leading-order chiral EFT that is the laboratory for 
our analysis. In Section III we examine decoupling for 



the deuteron and phase shifts starting from the LO EFT. 
The implications are discussed in Section|IV]and we sum- 
marize in Section IVl 



II. BACKGROUND 

A. Chiral EFT at Leading Order 

An implementation of chiral EFT at leading order 
(LO) in Weinberg's power counting was studied in 
Ref. [12 for a large range of momentum cutoffs Aeft 
to determine what renormalization was needed in each 
partial wave to achieve cutoff independence. The EFT 
NN potential at LO consists of one-pion exchange in all 
channels plus a regulated contact interaction in certain 
channels; it is meant to be iterated to all orders. The 
power counting proposed by Weinberg included contact 
interactions in only the S- waves at LO, but this was found 



by Nogga et al. to be inadequate for cutoff independence 
in partial waves that have an attractive tensor force (for 
high enough angular momentum there is no problem). 
Adding a single contact interaction in applicable chan- 
nels remedied this problem and allowed large cutoffs to 
be explored. 

A side effect of the attractive tensor force was the ap- 
pearance with increasing AgpT of spurious deep bound 
states. These are not a problem in principle for the EFT, 
because they appear outside its domain of validity. That 
is, the EFT is expected to have incorrect ultraviolet (UV) 
behavior, whose impact on low-energy observables is cor- 
rected by the contact counterterms. In this case, UV 
refers to large negative as well as positive high energies. 
In practice the problems caused by spurious deep bound 
states for solving few-body equations led Nogga et al. to 
remove them by hand |12^. Our question here is: Can 
these deep bound states be decoupled by renormaliza- 
tion group flow equations? This question makes the LO 
chiral EFT Hamiltonian a good test case to study the 
robustness of the observations in Ref. PD] and to learn 
about the decoupling properties of the SRG as a function 
of Aeft- The present study is restricted to two-particle 
systems. The impact on three- and higher-particle sys- 
tems may be of even greater interest, but is deferred to 
future work. 
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FIG. 1. (color online) Dimensionless coupling constants for 
the LO chiral EFT as a function of the EFT cutoff Aeft- For 
the ^So channel, Cs = Ca/ is fit to the phase shift at Eiab = 
lOMeV. For the coupled ^Si-^Di channel, ct = ct/fl is fit 
to the '^Si phase shift at lOMeV. The constant Cp = Cp/ f!^ is 
fit to the ^Po phase shift at 50MeV. 

Wc follow Ref. [H] in defining the LO chiral EFT. For 
all our calculations, we regulated the interaction with 



/„(fc , k') = e^'*-' MEFT)^"g-(fc7AEFT)^" ^ 



(6) 



using n = 4. Our notation differs from Ref. |12| only 
in the definition of dimensionless coupling constants (see 
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Fig. [T] caption). We focus on three channels, which are 
sufficient for the present discussion. The ^Sq channel 
is an example of a partial wave where the tensor force 
does not contribute and therefore has no spurious bound 
states. The coupled '^Si-'^Di channel and the ^Pq chan- 
nel provide representative examples of partial waves with 
an attractive tensor force; each develops spurious bound 
states with increasing Aeft in addition to the physical 
deuteron in the coupled channel. 
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FIG. 2. (color online) Bound-state energies as a function of 
Aeft for the ^Si-'^Di coupled channel (deuteron as a solid 
line, spurious states as dotted lines) and the "^Po channel (spu- 
rious states as dashed lines). 
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FIG. 3. (color online) '^So phase-shifts as a function of labo- 
ratory energy for several different Aeft. 

For the singlet and triplet S-waves we fit the cou- 
pling constants to the neutron-proton phase shifts at 
-Eiab = lOMeV, and for the "^Fq channel we fit to the 
-E'lab = 50MeV phase shift. The couplings in these chan- 
nels are shown in Fig. [l] and the bound state energies 
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FIG. 4. (color online) ^Si phase-shifts as a function of labo- 
ratory energy for several different Aeft. 
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FIG. 5. (color online) ^Po phase-shifts as a function of labo- 
ratory energy for several different Aeft. 



in Fig. [2j each as functions of Aeft- We see that the 
running of the coupling constant in the ^Sg channel is 
quite mild throughout the range of Aeft- In contrast, 
the running in each of the other channels displays a limit- 
cycle pattern. With each cycle to -|-c» a new bound state 
emerges, at energies bound by two hundred MeV or much 
more, which place them outside the low-energy domain 
of the EFT. 

In Figs. [3] |4j and [5] we plot the phase shifts for repre- 
sentative values of Aeft up to a lab energy of 100 MeV. 
Because the EFT is only evaluated at leading order, we 
should not be surprised at large deviations compared to 
experiment (here represented by the Nijmegen partial 
wave analysis). However, our interest is not in estab- 
lishing how well experiment is reproduced, but in how 
close the predictions of the various Hamiltonians are to 
each other as a function of interaction energy. This will 
be relevant when we examine universality below. 
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FIG. 6. (color online) Contour plot of V\{k, k') in the ^So channel for xeft with a cutoff of Aeft = 8.0 fm ^ for Trd (top) and 
Hd (bottom) SRG evolution. 



III. RESULTS 

In this section, we evolve LO chiral EFT potentials at 
representative cutoffs using flow equations with genera- 
tors based on both Ti-ei and Hd- We start with the ^Sq 
channel, for which there are no spurious bound states. As 
illustrated in Fig. |6]for a cutoff of Aeft = 8.0 fm~^, the 
evolved potentials for T^.^i and Hd are practically indis- 
tinguishable. This is true at any cutoff for this channel. 
More generally we have found it to be true for any chan- 
nel for which there are no spurious bound states. 

A very different story as a function of Aeft is found 
in the coupled ^Si-^Di channel, as seen in Figs. [Tj [sj 
and [9j which compare the Tj-oi and Hd evolutions for 
Aeft = 4.0, 9.0, and 20.0 fm~^, respectively. The initial 
Aeft ~ 4.0 fm^^ interaction supports only the deuteron 
bound state. Without any spurious states to decouple, 
the evolutions generated by T^ei and Hd are nearly iden- 
tical to each other, as with ^Sq. However, in Fig. [sjthere 
is now a spurious state to decouple and so the SRG evo- 
lution for Tici and Hd are substantially different. 

Both the Tj-oi and Hd evolutions steadily drive the ma- 
trix toward a diagonal form. In the process, the low 
momentum part of the matrix is driven towards large 
negative values. With T^c\ this is unabated, while in 
the Hd evolution it reverses momentarily and deposits 
an isolated negative eigenvalue on the diagonal, which 
corresponds to the complete decoupling of this spurious 
state (201 . The low- momentum potential in the last two 
panels of Fig. [s] for Tid (bottom) is qualitatively differ- 
ent from the corresponding panels for Aeft = 4.0 fm^^ 
in Fig. [7| while the corresponding Hd panels are similar 
for momenta below the spurious state matrix element. 

Figure. [9] shows what happens when a second spuri- 
ous state is supported after increasing the EFT cutoff to 
Aeft = 20.0 fm"\ A gain Troi and Hd differ significantly 
but the latter is similar to Hd in the other figures. If 
we expanded the momentum scales, we would find that 



the deepest state is first decoupled, with the diagonal el- 
ement found at much larger k, and then the shallower 
spurious state is decoupled at smaller A. 

For the uncoupled '^Po channel, the tensor is also at- 
tractive and increasing the cutoff of the EFT leads to a 
spurious bound state for Aeft between 3 and 4 fm""'^. In 
Figs. [TOj [TlJ and [12] we trace the evolution for two cut- 
offs above the threshold for a single bound state (Aeft = 
4 fm^^ and Aeft = 8 fm~^), and for a third above the 
threshold for two bound states (Aeft = 20 fm~^). Re- 
sults for low cutoffs with no bound states are not shown, 
but have the same pattern as for ^So. With a spuri- 
ous state to decouple, the SRG evolutions for Tj-ci and 
Hd at the lowest Aeft are again substantially differ- 
ent, just as for the '^Si-^Di channel. We see particu- 
larly non-universal evolution of the low-momentum part 
of the Tj-ei potentials while the Hd potentials are similar 
(more quantitative comparisons are made below). Note 
that the diagonalization of the spurious bound states do 
not occur at the same diagonal momentum. 

The quantitative evolution can be characterized by fo- 
cusing on the diagonal and fully off-diagonal matrix el- 



ements separately. In Figs. 13 and 14 the diagonal and 
fully off-diagonal matrix elements of the interactions in 
the '^Si are plotted for Hd and T^-c\, respectively, as rep- 
resentative examples. The decoupling of high and low 
momentum degrees of freedom is a prevalent feature in 
the Wegner evolution whether there is a deep bound state 
or not. The result of this decoupling is that the low mo- 
mentum portion of the matrix becomes universal, with 
a collapse of the original potentials to almost the same 
low-fc dependence for sufficiently low A. In contrast, with 
the Ti-oi evolution the low-momentum diagonal matrix el- 
ements are always quite distinct. The plots showing the 
off-diagonal edge manifest the distinction most clearly 
because the non-universal off-diagonal elements above A 
are driven to zero. Similar behavior is observed in other 
channels. The degree of collapse with the Wegner evolu- 
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FIG. 7. (color online) Contour plot of V\{k, k') in the "*Si channel for Xeft with a cutoff of Aeft = 4.0 fm ^ for Tj-d (top) and 
Hd (bottom) SRG evolution. 
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FIG. 8. (color online) Contour plot of Vx{k, k') in the ^Si channel for xeft with a cutoff of Aeft = 9.0 fm ^ for Tj-d (top) and 
Hd (bottom) SRG evolution. 
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FIG. 9. (color online) Contour plot of Vx{k,k') in the ^Si channel for Xeft with a cutoff of Aeft ~ 20.0fm ^ for Troi (top) 
and Hd (bottom) SRG evolution. 



6 



0~ 

ri 

E 





k (fm-i) 

,0 1 2 



k (1m-i) 

1 2 



k (fm-i) 
1 2 



k (fm-i) 

1 2 



k (fm-i) 

1 2 






■D.O{fm) 



^20.0 



A = 3.0 fm- 



A = 2.8fm- 



A = 2.0 fm- 



A = 1 .2 fm- 



A = 0.8 fm- 



FIG. 10. (color online) Contour plot of V\{k,k') in the "'Pq channel for xeft with a cutoff of Aept ~ 4.0 fm ^ for Troi (top) 
and Hd (bottom) SRG evolution. 
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FIG. 11. (color online) Contour plot of V\{k,k') in the ^Pq channel for xeft with a cutoff of Aeft ~ 8.0 fm ^ for r^ci (top) 
and Hd (bottom) SRG evolution. 
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FIG. 12. (color online) Contour plot of Vx{k,k') in the "^Pq channel for xeft with a cutoff of Aeft ~ 20.0 fm ^ for T^^i (top) 
and Hd (bottom) SRG evolution. 
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FIG. 13. (color online) The diagonal matrix elements Vs{k, k) and fully off-diagonal matrix elements Vs{k, 0) for the ^Si channel 
are plotted against momentum k after evolving using the Wegner generator to a range of final flow parameters A = s~^^*. For 
each A, results for flve cutoff values Aeft are shown (in some cases they are beyond the limits of the plot). 




k [fm '] k [fm '] k [fm '] k [fm '] k [fm '] 

FIG. 14. (color online) The diagonal matrix elements Vs{k, k) and fully off-diagonal matrix elements Vs{k, 0) for the '^Si channel 
are plotted against momentum k after evolving using the Trd-based generator to a range of final flow parameters A = s~^^^. 
For each A, results for flve cutoff values Aeft are shown (in some cases they are beyond the limits of the plot). 



tion is correlated with the level of agreement of the phase 
shifts in that channel. 



Figure 15 shows similarly striking consequences of the 
choice of generator for the momentum distributions of 
the bound states in the '^Si-'^Di channel, with Aeft = 
20 fm~^ and an evolution to A — 1.2 fm~^ chosen as a 
representative case. The initial deuteron momentum dis- 
tribution has a nearly exponential tail while the spurious 
state is essentially flat in the region of momenta plotted 
(reflecting the confinement of the coordinate-space wave- 
function). Their fates after evolution with Troi and Hd 
are starkly different. In the former case (left panel), the 
spurious state has had its strength pushed to low mo- 
mentum while the deuteron becomes broadly distributed 
and qualitatively changed at low momentum. In the lat- 
ter case (right panel), the spurious state is confined to 
a small region in momentum near k — 1.6 fm~^ while 
the deuteron momentum distribution at low momentum 



is essentially unchanged. Similar results are found with 
other choices of Aeft- 

The phase shifts for a given Aeft are unchanged as a 
result of a unitary SRG evolution. However, we can test 
the degree of decoupling in different cases by cutting the 
evolved potential (i.e., setting it to zero for all k, k' above 
a cut momentum) and looking at t he degrada tion in the 
phase shifts. This is shown in Figs. 16 17 18 for the '^Si 



phases shifts evolved from Aeft = 4, 9, and 20 fm , re- 
spectively. The cut is implemented by multiplying each 
matrix element by the factor e 



-{k/Ay 



Re- 



sults for m = 2, 4, 8 are similar; we show m = 4 only. In 
general, successful decoupling means that the phase shifts 
calculated with a potential evolved to A are unchanged 
at low energy if A < A (of course they will deviate for 
energies above the cut). 

If there are no spurious bound states, decoupling works 
for either T^^i or Hd evolution, as seen in Fig. 16 In 
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FIG. 15. (color online) Momentum probability densities in the coupled ^Si-^Di channel with Aeft = 20 fm^^ for the deuteron 
(d) and shallowest spurious state (sp). On the left are the initial and evolved densities using the Troi generator evolved to 
A = 1.2 fm~^ while on the right are the same densities but evolved using the Hd generator. 



particular, the phase shifts are essentially unchanged for 
A < 2.0 fm"^ after a 2.5 fm~^ cut and agree at low en- 
ergy after a 1.5 fm^^ cut for sufficiently low A. The 
pattern of decoupling is similar for i?d-evolved poten- 
tials when there are spurious bound states, except that 
the onset of decoupling is delayed until A is sufhciently 
below where the bound state is deposited on the diag- 
onal. Thus, in the top panels of Figs. [17] and [18) de- 
coupling is not achieved until A = 1.2 fm~^, compared 



to A = 2.0 fm"' in Fig. 16 for the 2.5 fm"' cut. This 



may slow convergence in few- or many-body calculations; 
this remains to be investigated. The decoupling patterns 
for phases shifts calculated from Trci-evolved potentials 
with spurious bound states show serious distortions at 
the lower cut. sufficiently low the phase shifts are again 
reproduced at low energy. 



IV. DISCUSSION 

We have documented several different outcomes from 
evolution with generators using and T^ei- For T^■c\^ 
the different classes of low-momentum Hamiltonians de- 
pend on the number of spurious states, while for 
there is a single universal class for the low-energy part, 
after sufficient evolution. There is corruption of the 
low-momentum probability distribution for Tid and the 
physical low-lying bound state is still coupled and non- 
universal. In contrast, with Hd the spurious state com- 
pletely decouples and the deuteron wave function is un- 
changed at low momentum. Note that despite these 
unhappy features for Tj-^i, the results for all observ- 
ables are formally the same because we still have a uni- 



tary transformation acting on a renormalized Hamilto- 
nian. In practice, however, numerical computation may 
lose all precision for observables influenced by the near- 
diagonalized part and our ability to truncate a basis ex- 
pansion will be compromised. 

Although it is clear that decoupling increases with de- 
creasing A when using Hd, it may not be so useful in 
practice. The binding energies of the spurious states seen 
here are very large, yet the decoupling momenta for the 
least bound ones are still comparable or smaller than 
values where SRG evolution is usually halted because of 
concerns about the growth of many-body forces. Will 
the remaining distortions spoil the advantages? For ex- 
ample, are variational calculations still effective and is 
convergence in basis expansions still accelerated? An- 
swering these questions will require extending our tests 
to A > 3. Also, the question of where a bound state is 
placed on the diagonal in momentum representation is 
still open. 

In Sect, jn] we invoked the claim that the deep bound 
states are outside the range of the EFT and as a result 
should not disturb the low-energy EFT predictions. Our 
results provide an explicit test of this principle. In partic- 
ular, the SRG unitary flow does not change the S-matrix 
and yet with an iJ^-based generator we are able to create 
a decoupled version of the Hamiltonian. Because it be- 
comes decoupled at a finite momentum scale A and has 
the same universal form in the low-energy region as po- 
tentials with no spurious bound states, we have achieved 
the desired demonstration. 

An important open question is the impact on SRG evo- 
lution of bound states for A > 2. For example, at what 
A do we expect to have problems evolving with a Tj-ci- 
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FIG. 16. (color online) Phase shifts in the coupled '^Si-^Di channel using the Aeft = 4.0 fm ^ potential, testing Ha (top) and 
Trei (bottom) SRG decoupling by first cutting off the potential. 
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FIG. 17. (color online) Phase shifts in the coupled ^Si-^Di channel using the Aeft = 9.0 fm ^ potential, testing Hd (top) and 
Trci (bottom) SRG decoupling by first cutting off the potential. 
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based generator for the triton, ^He, or nuclear matter? 
The choice between Tid and Hd for building the SRG 
generator is just one of many possible useful choices to 
control the flow of the Hamiltonian. What other gen- 
erators achieve what Hd does but allow more freedom? 
Work is in progress to address these issues. 



V. SUMMARY 

In this paper, we have applied SRG evolution flow 
equations that have been uniformly successful with com- 
monly used initial NN interactions to a new class of 
initial potentials. These potentials are from renormal- 
ized leading-order chiral EFT and have much higher 
momentum-space cutoffs which, combined with singular 
pion exchange, lead to spurious bound states emerging 
in channels where there is a non-zero, attractive tensor 
force. We have tested flow properties such as decoupling 
and the approach to universal form. 

While the choice of in the SRG generator has been 
successful for nuclear applications to date, this has been 
the fortuitous result of shallow bound states that are 
not disturbed for the A's considered. The presence of 
deep bound states introduces problems first emphasized 
by Glazek and Perry [20]. The generalization to Hd re- 
stores the good features such as decoupling and flow to 
a universal form as A decreases. 

An important moral is that even a comparatively sub- 
tle change in the generator can have a substantial ef- 



fect on the flow of the Hamiltonian. While in principle 
observables are always unchanged because we are mak- 
ing unitary transformations, in practice there are always 
truncations that make the transformations only approx- 
imately unitary. It is critical to minimize the impact of 
these approximations as we extend calculations to more 
particles. At the same time, the degree of perturbative- 
ness and convergence properties of the Hamiltonian can 
also be significantly affected. Future work includes look- 
ing at few-body systems and using the SRG in further 
investigations of EFT. 
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